Nonanalyticity of the optimized effective potential with finite basis sets 
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We show that the finite-basis optimized eff'ective potential (OEP) equations exhibit previously 
unknown singular behavior. Imposing continuity, we derive new well-behaved finite-basis-set OEP 
equations that determine OEP for any orbital and any large enough potential basis sets and which 
adopt an analytic solution via matrix-inversion. 
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I. INTRODUCTION 

In the last couple of decades the optimized effective 
potential (OEP) theory [l|-lJl appeared to offer a very 
promising route for iinproved accuracy in density func- 
tional theory (DFT) 0-0]. With OEP, not only is the 
exchange potential determined exactly but also there is 
hope that correlation will eventually be approximated ac- 
curately, via an implicit density functional |8l-lld|. Fur- 
thermore, using OEP and imposing physical constraints 
on the potential, it is possible to improve the performance 
of traditional approximations, such as the local density 
approximation [ll| . 

Recently, it was discovered that finite basis implemen- 
tations of OEP arc marred by mathematical problems 
[13 | . Several attempts have been made to overcome these 
issues J13l - [2G | but with limited success so far. As a result 
and despite promise, interest in OEP has diminished. 

The OEP is determined by a Frcdholm integral equa- 
tion of the first kind. 



dv'x,{v,v')v{v')^K{v). 



(1) 



In this work, the OEP u(r) represents the sum of the 
Hartree and exchange-correlation potentials in the effec- 
tive single-particle Hamiltonian, 



hy{r) = — — + Vc„(r) + v{r) , 



(2) 



where Won(r) is the electron- nuclear attractive potential. 
The single-particle eigenfunctions of h^ and their energies 
build the density-density response function in ([IJ, 

(r r') -2 V '^"-''^'^''^"''''^'^-''^''-"'-'^''''^"-' '-''''' (3) 
^- — ^ ^. ■ — ^. _ 



Indexes i and a run respectively over occupied and un- 
occupied orbitals in the OEP Slater determinant. &„(r) 
in (m has units of density. For concreteness, we omit 
correlation and focus on exchange OEP (x-OEP), where 

6,(r) = 2 Y. ^"'^""^^'^^ 0.,(r) ^,,Jr) , (4) 



where Ji, (r) is the direct Coulomb (or Hartree) local po- 
tential operator and ICy is the Coulomb exchange nonlo- 
cal operator. We use the shorthand {i\-\a) for the matrix 
element ((^„,i| • \(j>v,a)- 

Hirata et al. [lj| proved that, for a complete or- 
bital basis set, including continuum states, the products 
<f>T}^i{v) 01,. a (r) form a complete set bar a constant. Hence, 
the response function ^ is defined over the whole vec- 
tor space of functions orthogonal to the constant func- 
tion. Assuming that the response function is invertible 
in its space of definition, i.e. that the null-space of Xv 
is only the constant function, we conclude that OEP is 
fully determined up to a constant. 

On the other hand, with a finite orbital basis set the 
straightforward search for OEP may yield an infinity of 
solutions [12|: from Eqs. (II])-(|11) follows that the potential 
is undetermined in the null-space of the response func- 
tion. The latter space contains only the constant func- 
tion in the case of full Xv- With a finite orbital basis 
set the infinite sum over virtual orbitals in Eq. ([3|) is re- 
stricted to a sum over a subset of the virtual orbitals, in 
particular only over those obtained by the orbital basis 
functions. The truncated response function x° obtained 
in this way has an infinite-dimensional null-space. This 
results in indeterminacy of the potential along any of its 
components that lie in the null-space of x% ■ 

Interestingly, several approximations of the finite-basis 

OEP approach |2l| - [25j . which invariably employ the 

Unsold approximation [26|, determine the approximate 
OEP fully. The Unsold approximation amounts to a com- 
mon energy denominator approximation for the static or- 
bital Green's function, together with the closure or com- 
pleteness relation. The reason is that with the closure 
relation the orbital basis set becomes effectively complete 
and the null-space of Xv reduces once more to the con- 
stant function. Consequently it is no longer possible for 
an auxiliary basis function to have a component in the 
null-space of Xv ■ 

Returning to the OEP, the remedy appears readily, at 
least in principle. Consider a finite or even a complete 
auxiliary basis set {Cn(r)} for the expansion of the po- 
tential and take the matrix elements of the truncated 



response function Xv with the auxihary basis functions: 

A,r.^-lJJdrdr'Ur)xl{r,r')Ur')- (5) ^^ 

This matrix gives the projection of xS i^^ the space 
spanned by the auxihary basis functions. Its diagonahza- 
tion yields eigenfunctions with singular and with nonsin- 
gular eigenvalues, which provide a convenient orthonor- 
mal basis in the auxiliary space. The null eigenfunc- 
tions of Akn are auxiliary functions that belong in the 
null-space of Xv ^^'^ cannot be used for the expansion of 
the potential, since the potential is undetermined along 
them. The expansion of the potential along the remain- 
ing nonsingular eigenfunctions is well defined. 

Unfortunately, such a singular value decomposition 
(SVD) and truncation of the null eigenvectors is in gen- 
eral ambiguous because the singular eigenvalues of the 
matrix Ai^n are not always separated unambiguously 
from the nonzero ones, as can be seen in Fig. [1] for the 
correlation-consistent polarized-valence triple-zeta (cc- 
pVTZ) basis sets. Then, the inversion of Akn is not 
straightforward. 
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Figure 2: Finite-basis exchange OEPs v^ for the Ne atom 
with the same cc-pVDZ orbital basis set and four different 
auxiliary basis sets: uncontracted cc-pVDZ, uncontracted cc- 
pVTZ, uncontracted cc-pVQZ, and uncontracted Partridge. 
The four potentials converge with an expanding auxiliary ba- 
sis set to a potential that depends on the orbital basis only. 
The full numerical result from Ref. [l^ is shown with a dotted 
line. 
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Figure 1: The eigenvalues of matrix A for the Ne atom. Cir- 
cles, squares, diamonds, and triangles correspond to the cc- 
pVDZ orbital basis set and uncontracted cc-pVDZ, cc-pVTZ, 
cc-pVQZ and Partridge auxiliary basis sets respectively. Stars 
correspond to cc-pVTZ orbital and uncontracted cc-pVTZ 
auxiliary basis sets. 

Nevertheless, sometimes, the singular eigenvalues of 
Akn can be unambiguously identified a priori, with a 
clear gap of many orders of magnitude separating them 
from the rest. Such a case is shown in Fig. [T] with the 
correlation-consistent, polarized-valence double-zeta (cc- 
pVDZ) orbital basis set and several auxiliary basis sets. 
In this case, the resulting potential after truncation of the 
singular eigenvectors and inversion of Akn is expected to 
be unique. However, it is a mystery why in these cases 
the resulting potential, which is mathematically unique 
and well defined, looks unphysical with strong oscillations 
appearing near the nuclei; in the case of atoms, these os- 



cillations make the potential look very different from the 
exact x-OEP [l^, as can be seen in Fig. [21 

In Fig. [2] we show the x-OEPs for the same orbital 
basis set cc-pVDZ and four auxiliary basis sets cc-pVDZ, 
cc-pVTZ, cc-pVQZ and Partridge (all uncontracted). As 
expected, since the truncation of the null eigenvectors 
of Akn is clear and straightforward, the four potentials 
in Fig. [2] are close to each other and converge with an 
expanding auxiliary basis set. Thus, Fig. [2] provides 
a numerical demonstration of our argument that when 
the SVD of matrix Akn is unambiguous, the resulting 
potentials from finite-basis OEP theory are mathemat- 
ically well defined. Of course, it is surprising that the 
plots converge to a potential that looks unphysical and 
is significantly different from the full numerical x-OEP. 

So far in the literature, the underlying reason for this 
anomaly has remained elusive and it is usually confused 
with the ill-posedness of the inversion of Akm the two 
problems appearing as one. For example, it is often writ- 
ten that similar wild oscillations of the potential near 
the nuclei appear when the orbital and auxiliary basis 
sets are not balanced. The lack of precise definition of 
balanced basis sets notwithstanding, we have observed in 
Figs. [1] and [2] that the anomaly is present after the trun- 
cation of auxiliary functions in the null-space of Xv > even 
in cases where the SVD of Akn is clear and unambiguous. 

In this work, we distinguish between the cause for 
the anomalous behavior of the otherwise mathemati- 
cally well-defined finite-basis OEP, and the general ill- 
posedness of the choice of SVD cutoff for the truncation 
and inversion of Akn- We shall first focus on the for- 
mer clear-cut problem which we shall resolve. Then, we 
shall discuss briefiy the more common and complicated 



case, where the anomaly appears entangled with the ill- 
posedness of the inversion of the matrix A^^ of the re- 
sponse function. Wc shall investigate this problem at 
length in a future publication. However, preliminary re- 
sults presented in the last section indicate that our treat- 
ment of the nonanalyticity helps to sort out this problem 
as well. 



II. ANALYSIS 

We shall base our analysis on the solution of the inter- 
mediate x-OEP equation, 



X.U '^oep — ^u : 



(6) 



with a fixed potential u in ft,„ [Eq. ([2])], i.e., u is not 
updated to self-consistency; we also omit the explicit de- 
pendence on u. For simplicity, we shall analyze the con- 
sequence of a finite orbital basis set in the special case 
where the finite orbital set is composed of the occupied 
orbitals (pi and of a subset of the (mostly lower lying) vir- 
tual orbitals (pa of h. The virtual orbitals of h [Eq. ^] 
that are outside the orbital basis form the complement 
basis; they will be denoted by (pa (or simply a) and their 
energies will be denoted by la- 
in the first part of our analysis we consider that the po- 
tential is described on a grid with arbitrary precision, or 
equivalently that the potential is expanded in a complete 
auxiliary basis set. 

To study the effect of the finite orbital basis we split 
X and h in Eqs. (H])-®: 



X(r,r') = Xo(r,r') +x(i",r'), 



h{v) = bo{v)+h{v), 



(7) 



(8) 



where xo ^-nd &o a-rc given by Eqs. (jH]) and Q but the 
sums over virtual states are restricted in the orbital finite 
basis. The remainder functions (with a tilde) are defined 

by 



X(r,r') - 2 5] 
6(r) = 2^ 



0,(r)^,(r)0,(r')0,(r') 



{h\J -K.\i 



(9) 



^i - Eq 



Cp^{v)cPa{v). (10) 



The summation over virtual orbitals is in the complement 
basis. 

Let us denote by v^ the potential which satisfies the 
equation, for A > (we drop the subscript oep), 



(xo + Ax)w^ = fco + A6. 



(11) 



The exact x-OEP is obtained for A = 1. In a finite-basis- 
set implementation the unknown x is always omitted. 



which amounts to setting A = 0. Equation (TTTI) reduces 
to 



Xow 



with 



Xo 



'b. 



(12) 



(13) 



in place of w^, representing the finite-basis x-OEP solu- 
tion. 

However, we point out that for finite A > 0, the re- 
sponse function (xq + A x) in (|lip is invertible and v^ is 
fully determined up to a constant. On the other hand, 
for A = the response function reduces to xo: which has 
an infinite-dimensional null-space where v^ is undeter- 
mined. Hence, the complete omission of x is a singular 
operation and the following theorem holds: 

Theorem 1. The solution v^ of the OEP equation pTJ) 
is not a continuous function of A at A = 0. 



''^v\. 



(14) 



where v stands for limA-yo "^ ■ 

To expand on the proof of the theorem, we consider 
the eigenfunctions c" (r) of xo with nonzero eigenvalues 
and an orthonormal basis c'^(r) in the null-space of xo^ 
(i.e. g- = 0) 



ydr'xo(r,r')c"(r') = -2.g"c°(r) 



(15) 



The index a enumerates the nonsingular eigenfunctions 
and v, ^ run over the null eigenfunctions. The functions 
c''(r) cannot be chosen to be eigenfunctions of x since the 
latter has a projection in the nonsingular space of xo- 

We shall use the complete set of eigenfunctions of xo 
as a special basis in which we expand the potential. 

From Eqs. ^^ and ([T3|), the potential v°{r) is ex- 
pressed in terms of the nonsingular eignfunctions c"(r): 



."(r)=^ 



6„c°(r) 



where ba is the overlap: 

6„ = -iydrc"(r)foo(r) 
We also define the overlap: 



K^-- I drc%Y)h{v). 



(16) 



(17) 



(18) 



We expand in Taylor series v^, around the value w", 
with A > K > 0: 



w-^ ^ v'' + {\-k)v'^ + 



> , i^~^)\.>( 



V' + 



(19) 



We take the limit k — )► 0. To allow for the possibility 

^K-)-o ^ 7;0^ tfjQ write: 



v^^v^^^+Xv' +'—v" + .., 



A2 
2 



(20) 



To determine v^^°, we substitute (HO]) in dlT]), use (IT2l) . 
and keep up to first order in A: 



Xqv 



A^O 



+ X{Xov' + Xv^^'')^Xov'' + ^b. (21) 



From the zero-order equation, xo v^^^ = xo i'^ , we ob- 
tain that V = v^^^ — ti° is a null vector of Xo and can be 
expanded in the null eigenfunctions of xo j 



i;(r)=^i;^c^(r). 



(22) 



(23) 



The potential ii is a measure of the discontinuity. 
The linear term in Eq. ([21]) gives: 

Xov' + xv" + xv^b. 

Multiplying by c'^(r), integrating over r, and using 

ydrc''(r)xo(r,r') = 0, (24) 



we obtain, 



6a 



/ ^ Xvn '^'/J. — bi, 2_^ Xi^a 



(25) 



where. 



X.^ = -2 // drdr'c''{r)xir,r')c^{r'), (26) 



X.^ 



dr<ir'c''(r)x(r,r')c"(r'). (27) 



The null-space of xo (where c'^ and c'' lie) is a proper sub- 
set of the nonsingular space of x and hence the matrix 
X^p, is well-defined and invertible. Also, the nonsingu- 
lar space of X overlaps with the nonsingular space of xo 
(where c" lie) and the matrix Xz^q does not vanish iden- 
tically. We finally have: 

^t^=Yl ^'t^^ ^"^ " 51 ^M^^ ^-^^ -^ ■ (28) 

In general, the right-hand side and hence v are not ex- 
pected to vanish. This completes the proof that v^ is 
discontinuous at A = 0. QED 

Equation ([28]) determines v in the null-space of xo- In 
order to calculate v from Eq. ([28]) we need some knowl- 
edge of X- In the next section we shall use the Unsold 
approximation to approximate x- 

To conclude this section, it is important to note that 
the nonanalyticity of OEP at A = is the result of the 
truncation of the orbital basis set. So far we have taken 
that the potential is either represented on a grid with 
arbitrary precision, or that it is expanded in terms of 
the complete set of eigenfunctions of xo- Therefore, the 
finiteness of the auxiliary basis set (in practice) has not 
played a role. In the following, we shall expand the po- 
tential in an arbitrary but complete auxiliary basis set to 
complete the derivation of the new finite-basis-set OEP 
equations in which the potential is not discontinuous at 
A = 0. 



Arbitrary complete auxiliary basis set. 

So far we have expressed the potential in the basis 
of the eigenfunctions {c"(r)} of the truncated density- 
density response function xo(r,r')i complemented for 
completeness by an orthonormal basis {c'^(r)} in the null- 
space of xo(r,r'). Although {cP(r)}, (p = a, i^), has been 
a natural basis to base our analysis on, it is not a practi- 
cal basis for calculations and in this section we shall make 
a change of basis in the auxiliary space and expand the 
potential in an arbitrary but complete auxiliary basis set 
{^„(r)}. Of course, in practice, the auxiliary basis set is 
never complete; however, the finiteness of the auxiliary 
basis set is not expected to introduce any further nonana- 
lyticity in the potential. In addition, the new finite-basis 
OEP equations, derived in this section, give meaningful 
results for a finite auxiliary basis as long as the latter is 
large enough to overlap with the null-space of xo- 

In the following, we use a complete auxiliary basis set 
{£,n} to expand the eigenfunctions of the response func- 
tion. 



cP(r) = 5^cP^„(r), p = a, 



(29) 



Equation (P5)) defines the transformation between the 
auxihary bases: {c^(r)} -> {Cn(r)}. 

Substituting Eq. (P5)) in the eigenvalue equation (fT5|) 
and using Q, we obtain that the coefficients satisfy the 
generalized eigenvalue equations 



n n 

n 

where A^n in Eq. (fS)} is given explicitly by 

{^\ik\a){a\^n\^) 



Akn = Y^ ■ 



(30) 

(31) 

(32) 



^a ^i 



To distinguish from the previous basis of the eigenfunc- 
tions of xoj we use capital letters for the matrix elements 
of the response function with the auxiliary basis {^ri(r)} 
and for the overlap of the auxiliary basis functions with 
bo{r) and 6(r), 

Akn = -^ ffdrdr'^k{r)xir,r')Ur'), (33) 



Bk = -2 ldr^k{r)bo{r), 



Bk = --JdrUr)bir). 

From Q Bk becomes explicitly 

mk\a){a\J-IC\i) 



B.-Y.- 



ea - £» 



(34) 
(35) 

(36) 



Obviously, we cannot have x s-nd b exactly. We approx- 
imate [2g| the energy differences in the denominators of 
X and & by a constant, A ~ e„_a — Cu.i and then use the 
closure relation. The matrix elements become: 

Akn = E^*l^^-^"l*) - J2^i\^'^\^^(Mn\l) (37) 

Bk = E(*l^''(^ - ^)l*) - J2i'\^k\3){j\J ~ m (38) 

In Eqs. ([57]) and ([55]) we have ignored a term in Afe„ and 
a term in _Bfe, whose contribution in the OEP equation 
below, in the limit A — > 0, vanishes smoothly. 
Next, we expand the potential in the new basis, 



nr) = E""^"W- 

n 

In matrix form the OEP equation (|lip becomes 
Y^{Akn + A ifc„) vl = Bk + XBk, 



(39) 



(40) 



With definitions (|37l38p . A in Eq. ^ stands for A/A. 
In the limits A — s- and A ^> cxd our results are indepen- 
dent of A. 

For fixed A the solution of Eq. pOjl is 

(41) 



v::^J2(^ + ^^^nliBk + xB,) 



For A — )- oo, we observe that v°° satisfies the effective 
local potential (ELF) equations |24|.|25|: 



Y^^nlBk. 



(42) 



In the literature, the ELP potential v°° is also 
known as common-energy-denominator ajDproximation 
(CEDA)[23, or local Hartree-Fock (LHF)|22J. 

The same analysis, Eqs. ((20)) - ([28|) . which led us to con- 
clude that the potential v^ in Eq. pT|) is discontinuous 
applies to Eqs. (UHl) and (HJ). For small A, the potential 



tends to the limit 



,A^0 



v° + V. (43) 

The two components v^ and v are given by Eqs. p6p and 
We expand v^ and v in the auxiliary basis set, 



v\v) = Y^vlU^), v{r)=Y.^^Ur) 



(44) 



The change of the auxiliary basis {d'{r)} — > {^n(r)} is 
given by Eq. p9)) . It is straightforward to obtain 



a. 


(45) 


jj. 


(46) 


6„ = > Ic-^Bk, 

k 


(47) 


K = y]ciBk. 


(48) 



Finally, Eqs. dTe]) and ^^ become 






(49) 



'^n / , '-n /^fiv ^k -^k / ^ ^n A/^y Ai^a ^ 

l-iM,k oi,fj,,v ^ 

(50) 
In contrast to Eq. (pS)) . the matrices x^i/ and x^a above 
are not given exactly but within the Unsold approxima- 
tion: 

Xl-iu = /_^ C^ Akn C„ , Xva = 2_^ ^k ^-kn C„ (51) 



Therefore, Eqs. P^ and ([50)1 determine u° exactly but v 
only approximately. 

The finite-basis x-OEP equations gS]), gH), (gg]), and 
([50]) are the main result of this work. In these equations, 
the discontinuity of the potential v^ as a function of A at 
A = is evident. Until now, the potential w° played the 
role of finite-basis x-OEP. It is obtained by truncating the 
singular eigenvectors of A and subsequently inverting A 
in the space of nonsingular eigenvectors. The component 
V of v^~^'^ determines x-OEP in parts of the auxiliary 
space where v^ is undetermined. For example, the first 
term on the right-hand side of Eq. ([50)) is the projection 
of v°° in the null-space of Akn- Also if the null-space of 
Akn artificially included all the eigenvectors of Akn , then 
from Eqs. dH]), (gH]), and ^5U^, we would have v° = and 
t>^->o =v = v°°. 

The total energy as a function of A is continuous, i.e., 
the total energy of v^~^'^ is the same as that of ?j°. 



III. NUMERICAL EXAMPLES 

The eigenvalue spectrum of Akn for the Nc atom, using 
different basis set combinations for orbital and auxiliary 
basis sets, is shown in Fig. [1] In our calculations, we use 
Hartrec-Fock (HE) orbitals in place oi(j)u^i, eu,i, 4>u.a and 
eu,a- A gap separates the zero from the nonzero eigenval- 
ues for the cc-pVDZ orbital basis set and four auxiliary 
basis sets. In Fig. ^ we show that the four different po- 
tentials v'^ converge with expanding auxiliary basis to a 
potential that is characteristic of the orbital basis only. 
In Fig. [21 we perform the same test for the newly defined 
complementary potential v. Again, we observe that the 
four different potentials v converge and the converged 
potential depends only on the orbital basis. 

In Fig. m (Ne atom, cc-pVDZ and uncontracted cc- 
pVDZ orbital and auxiliary basis sets), we investigate 
numerically the interpolation of the potential v^ between 
the limiting potentials v^^'^ and v°° by plotting the po- 



tential difference v 



A 



,A^0 



The smooth convergence of 



i; toward the two limits for small and for large A can 
be seen clearly. On the scale of the plot, the potential 
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Figure 3; The exchange part, v^ of the correction poten- 
tials V, from Eqs. (|43ll50p . corresponding to the potentials v^ 
shown in Fig. [2] for the Ne atom. 



is on top of v°° already at A = 10. For A 



10- 



the potential v is indistinguishable from v and their 
difference vanishes. 
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Figure 4: Differences of the exchange potentials v^ — i^x^" 
for various values of A for the Ne atom with cc-pVDZ and 
uncontracted cc-pVDZ for the orbital and auxiliary basis set 
respectively. The potentials v^ are calculated using Eq. (|41|l , 
and Wx^° is calculated using Eqs. (gSJ, dM (EH and ([50)1 . 



The red long-dashed line is the difference v'^ — v^ 
v^ is calculated with Eq. (|i^ . 



where 



The limiting potential v^^'^ appears in Fig. EJa) (cc- 
pVDZ and uncontracted cc-pVDZ for orbital and auxil- 
iary basis sets), together with the potentials v^''^^^^ and 
v°° , and for reference the full numerical grid result from 
Ref. [H The potential v^^° lies almost on top of v°°°°^ 
and is quite close to the exact potential. Thus, by adding 
V to v'^ , the anomalous oscillatory behaviour of u*^ was 
corrected. This effect is further demonstrated in the inset 
of Fig.[SJa) where we show the two constituent potentials 



11° and V and their sum v'^^°. For this particular basis 
set, the total energy of x-OEP, for A = and A — ?■ 0, is 
identical to HF total energy. 

For each pair of orbital and auxiliary basis sets the 
plots of the potentials w" and v^^'^ are shifted by the 
same constant, which equals the difference between the 
highest occupied eigenvalues in the HF and the v^^'^ 



OFF calculations, 
potentials v. 



No shift is applied to the correction 



Ill-posed inversion of Akn 

In general, the problem of how to separate the sin- 
gular {g'^) from the nonsingular (g") eigenvalues of Akn 
is considered not straightforward and ill-posed [20, [231 ■ 
We focus on this problem in a forthcoming publication 
[28| . where we follow Ref. [20| and employ Tikhonov's 
regularization theory [29[, after deriving an appropriate 
regularizing function for x-OEP, from the energy differ- 
ence Tu [v] in Ref. |iO- In the rest of the section we present 
a simpler and preliminary analysis. 

An example where the logarithms of the eigenvalues of 
Akn spread almost uniformly between singular and non- 
singular values is shown in Fig. [1] for the Ne atom for or- 
bital and auxiliary basis sets cc-pVTZ and uncontracted 
cc-pVTZ, respectively. In Fig. (TJ we observe that the 
logarithms of the eigenvalues fall broadly in two groups 
with different slopes for increasing eigenvalue index p. In 
the first group, with p < 42, the slope is not steep and 
the logarithms of the eigenvalues fall off rather slowly. 
In the second group, with p > 43 the logarithms of the 
eigenvalues fall off faster. 

In Fig.[6l we plot the difference log^g (ff^^^) — logio (ff^) 
of consecutive eigenvalues versus the eigenvalue index p. 
We see there is a sharp increase in this difference as the 
eigenvalue index crosses the value from one group to the 
other. The lowest eigenvalue in the first group is g°' = 
1.04 X 10-3. 

Next, to obtain more rigorously the cutoff eigenvalue 
for the null-space of Akn , we start with trial values of the 
cutoff which are too large (with correspondingly too large 
trial null-space). We reduce the cutoff value gradually, 
changing the division between the effective null-space of 
Akn and the effective nonsingular space of Akn , by trans- 
ferring, one by one, eigenvectors from the former to the 
latter. This is straightforward since there is a finite num- 
ber of discrete eigenvalues. At each step we calculate the 
potentials, v{r) and w{r), from Eqs. (|M1), (SSI), and ((50)) 
before and after the shift of the eigenvector respectively. 
We monitor the step change of the potential by calculat- 
ing the root-square difference. 



Av = \ dr [v{r) — 'w{r)Y 



(52) 



/jj, denotes integration up to a large radius R, to pre- 
vent divergence of the integral. When the eigenvalues 




Figure 5: Exchange potentials, for A — >■ oo, A — >■ 0, of the Ne atom using two different combinations for the orbital and 
auxiliary basis sets: (a) cc-pVDZ and cc-pVDZ uncontracted and (b) cc-pVTZ and uncontracted cc-pVTZ. The full numerical 
result from Ref. [I3 on a grid is shown as exact. In the insets, the two strongly oscillating components of ifx"*" of Eq. (|28p are 
shown. Potentials for A = 0.0001 and A = 0.001 are also shown in (a) and (b), respectively, v^ in (b) is obtained by transferring 
an extra eigenvector (with eigenvalue < 1.04^^) from the effective null-space of A^n to the effective nonsingular space. 
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Figure 6: The step change of potential Aw (circles) and the 
quantity logj^Q (ff^^^) ~ logio id'') (stars) vs the index p of the 
eigenvalues (in order of decreasing magnitude) of the matrix 
Akn for Ne with cc-pVTZ and uncontracted cc-pVTZ orbital 
and auxiliary basis sets. 



are degenerate, we keep or remove the degenerate eigen- 
vectors together in each step. In Fig. [51 we plot the 
step change of the potential Av versus the eigenvalue in- 
dex p (Ne atom, basis sets cc-pVTZ and uncontracted 
cc-pVTZ, R=2 a.u.). It is shown that the resulting x- 
OEPs from Eqs. (gS]), (01]), and ([501) change slowly and 
converge until a truly singular eigenvalue is reached. At 
the step when the latter is transferred from the null-space 
to the nonsingular space, the change of potential Av is 
abruptly much larger than previous steps. We obtain 
that the smallest nonsingular eigenvalue of A^n is again 
g°' = 1.04 X 10-3. 

Comparing the two data sets in Fig. [51 it is evident 
that the behavior of the step change in the potential, ver- 



sus eigenvalue index p correlates well with the difference 
of the logarithms of consecutive eigenvalues, a simpler 
procedure that does not require the calculation of the 
potential in order to determine the cutoff. 

The corresponding exchange potential v^^'^ is shown 
in Fig. [5{b) for cc-pVTZ and uncontracted cc-pVTZ or- 
bital and auxiliary basis sets. We also show v^ (ELP) 



and 



i.e., the potential for A equal to the smallest 



nonsingular eigenvalue. Again, 



,,A->0 



and 



,,0.001 



lie al- 



most on top of each other. The oscillatory potential w* 
in Fig. Elb) results from transferring an additional eigen- 
vector (with the largest eigenvalue) from the null to the 
nonsingular space of Akn ■ The corresponding change of 
potential is the first large step change Av in Fig. [5] 

In the inset, as in Fig. Elja), we show the comple- 
mentary (oscillatory) potentials v'^ and v and their well- 
behaved sum w'^-*'". The total energies for ii^^o and 
^0.001 g^j.g almost identical and are 3 x lO^"* hartrees 
higher than the HF energy. In the last section, the no- 
tation "A -^ 0" merely impHes use of Eqs. (fT6)) . ((28l and 
(P5)) to calculate the finite basis x-OEP. with the cutoff 
for the null-space of A determined as described. Strictly, 
A cannot take values lower than the small but nonzero 
eigenvalues in the null-space of A. 



IV. CONCLUSION 

We have unveiled nonanalytic behavior of the OEP 
when the density-density response function is truncated 
with a finite orbital basis set. We proposed to employ the 
limiting potential v^^^ , instead of v^, for the appropri- 
ate description of finite-basis OEP. By using the Unsold 
approximation, we derived amended, finite-basis, x-OEP 
equations (|^5)) . (|491 and ([50]) . that determine v^^^ com- 
pletely for any combination of orbital and (large enough) 



auxiliary basis sets. The Unsold approximation amounts 
to employing an effectively complete orbital basis set, i.e. 
one that includes continuum states. Our new finite-basis 
OEP equations do not address the separate problem of 
how to distinguish between the effectively nonsingular 



and the singular eigenvalues of the finite basis matrix of 
the response function A^n- It seems however that with 
the help of the new finite-basis OEP equations, this tech- 
nical problem may become easier to tackle. 
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